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ABSTRACT 

This research develops concepts, flow diagrams, and statistical meth¬ 
odology l'or an information system to estimate the magnitude and scope ol 
nuclear attacks. The system is designed to operate in the transattack 
and immediate postattack period when data on the attack can be expected 
to be incomplete and inaccurate. The underlying principle in estimating 
attack si m is that of statistical inference, which permits an estimate 
to tie made of the total attack 1 com inlormation on only a sample oi the 
attack. Heuristic decision rules are applied as needed to make the system 
operable. It is demonstrated how the target synthesis procedure thus de - 
veloped may be coupled to a variety of survival estimating techniques to 
yield survival estimates. The research also develops statistical method¬ 
ology l'or processing reports of nuclear detonations. 
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of OCD'y propram to develop appropriate information systems relevant to 
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tistical methodology described in Section VIII of the re]xirt was developed 
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Wer* j devised by Frederic A. Miercorl. Benjamin t. Suta and Peter U. 
Bjorklund assisted in development ol the statistical methodology. 
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1 I f!TRODUCTI ON 


The purpose* of this report is to support research direct eel at improv¬ 
ing existing survival estimating systems, jirnl to develop new systems by 
which civil defense can undertake a realistic postattack assessment of the 
extent and damge of enemy attack and the resulting surviving people and 
resources. Executive Order 10952, as amended, assigns civil defense re - 
sensibilities in this urea as follows: 

"...develop plans ana operate systems to undertake a nationwide 
postattack assessment of the nature and extent of the damage re¬ 
sulting from enemy attack and the suiviving resources. Including 
systems to monitor and report s|»eeiflc hazards resulting from the 
detonation or use of special weapons..." 

In this investigation, effort Is centered on postattack assessment at 
the national level during the trnnsnttaok and immediate postattack periods. 
Postattack damage assessment systems at subnational levels are also con¬ 
sidered. but detailed investigations of these systems are left for future 

research 

The problem of estimating damage resulting from nuclear attacks on 
the continental United States has been the subject of continuing research 
for over a decade. During this time, researchers have devised numerous 
systems based on the use of high speed computers for determining damage to 
population and resources under almost any type of nuclear attack. The 
data produced by such systems have generally proved useful for purposes of 
pre ittack planning. On the other hand, survival estimating systems pro¬ 
viding information fer operations in a postattack environment have been 
studied to a much loss degree. In fact, a substantial part of the con¬ 
cepts and methodology of the damage assessment planning systems has been 
carried over for use with postattack damage assessment systems that are 
supposed to meet the information requirement for civil defense operations. 
In recent years, the continuing study of the nature of the problems in the 
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postattack environment has made it increasingly evident that the planning 
system damage assessment methodology will not supply adequate information 
for postattack civil defense operations. This is largely due to the tact 
that real-time information on the magnitude and scope of the attack will 
lie needed to apply existing damage assessment methodology, but will un¬ 
fortunately not be available. 

The accuracy and reliability of data available in the transattack and 
immediate postattack period may often be low. Attack and damage reports 
may be inadequate or biased. Nevertheless, the decisions that must be 
made duiing this period demand the availability of the best possible as¬ 
sessment of the extent of enemy attack. Accordingly, this investigation 
develops and describes an attack estimating system that will provide timely 
Information needed for operational purposes and that will minimize the ef¬ 
fects of inadequate and inaccurate information. It is further shown how 
the attack estimating system can lie used with existing damage assessment 
me t hods. 
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11 SUMMARY 


This research develops and presents a statistical concept for making 
national survival estimates in the transattack and early postattack 
periods. Statistical methodologies are used to process reports of detona¬ 
tions, and the principle of statistical inference is used to estimate the 
magnitude and scope oi the attack. tsting systems for the preparation 
of national survival estimates in the transattack and early postattack 
periods arc based largely on proattack planning methods or on the tabula¬ 
tion of direct damage reports as thc\ are received. Adequate and accurate 
data of the kind needed to produce survival estimates with these techniques 
will not be available in the transattack and early postattack periods. 

The four principal at tivitios oi the concept developed by this re¬ 
search are (1) real time input processing and data preparation, (2) sta¬ 
tistical conversion of data to information, (3) inference oi attack size 
and target system, (1) estimation of effects on resources and population. 
Activities 1 and d above may have much in common with certain existing 
survival estimating system-. The main distinctions between this and other 
survival estimating system- are in activities 2 and 3. The concept 
developed by this research u-o statistical inference and heuristic 
decision rules to estimate the magnitude of the total attack. The method¬ 
ology thus devised can operate with incomplete and inaccurate information. 
Targets arc inferred by statistical inference after appropriate grouping 
of target candidates into target categories and vulnerability groups. 

Tlie methodology includes: the receiving and processing of burst 
information and direct damage reports for use in making survival estimates; 
basic concepts for organizing inputs from various sources into appropriate 
data files; and provision lor supplying measures of reporting performance 
to the system operator. Detonation reports and other data are converted 
statistically to information that is more directly applicable to survival 
estimating. Estimated actual groun' ! zeros and targets for the initial 
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target list are determined statistically. Areas where data are incomplete 
or possibly inaccurate arc identified, and status reports are requested to 
correct these deficiencies. An initial target system is developed that 
yields information on the target categories hit ami on the severity of the 
a 11 ac k. 

The next main step is inference of the full attack size and target 
system, This is accomplished with the aid of the initial target list, a 
list of target candidates, and a knowledge of target categories and vul¬ 
nerability groupings. At the completion of this step, an augmented Target 
list that in general will be expanded considerably from the initial target 
list will have been developed through statistical inference. Again, 
status reports will be requested to cover areas of doubtful or inadequate 
information. 

The final step of the concept is to estimate effects on resources 
and population by means of the synthesized target system. It is shown 
iiow a target system may be coupled with several survival estimating tech¬ 
niques to produce survival estimates. The survival estimating technique 
selected may well depend upon the requirements oj the system operator at 
a given point in time. A summary of the system and its principal activi¬ 
ties is given in Table 1. 

Section VIII presents details of statistical methodology for esti¬ 
mating actual ground zeros and targets relevant to a group of detonation 
reports. A basic tool is the use of confidence regions and confidence 
intervals to determine the actual number of bursts and estimate actual 
ground zeros associated with a particular group of detonation reports. 

The selection of target candidates for the initial target list is also 
accomplished statistically, using a least-squares methodology. Weapon 
yields and heights of burst lor a given group of detonation reports are 
derived mathematically. 

Appendix A provides more detailed (low charts for some parts of the 
conceptual system. The basic purpose of these flow charts is to document 
the ideas thal have been developed in the course of this research but 
that are not appropriate to the generalized flow charts of the mam 
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report. The flow charts of Appendix A will help provide the basis for 
developing the system to a point where computer programs for the system 
can be written. However, considerable additional work will be needed 
before these flow charts are adequate for that purpose. 

The feasibility of Inferring detonation records from direct damage 
reports is investigated briefly in Appendix B, which also gives a frame¬ 
work for such inference, using pattern recognition and hypothesis testing. 

This rosearch effort has developed methodology for estimating attack 
size in the transattack and early postattack periods, and has shown how 
this methodology can be combined with existing damage assessment methods 
to produce national survival estimates. More advanced survival estimating 
techniques, that would integrate attack size estimating and damage assess¬ 
ment more directly, are briefly explored. Survival estimating systems 
recommended for full development and implementation by the National Civil 
Defense Computer Facility are identified. 

It was not the intention of this study to develop computer programs 
ready to use. The development of such programs, based on methodologies 
presented herein, will require much additional effort. 





Ill FUNCTIONS AND ORGA.NI ZATION OF CIVIL DEFENSE 

The general organization of civil defense is indicated in the Federal 
Civil Defense Guide,* which describes functions in the national, state, 
and local governments. The objectives and responsibilities of civil 
defense are outlined in the National Plan for Emergency Preparedness.t 

The principal objectives of civil defense are twofold: (1) to protect 
life and property by providing the means and knowledge necessary to mini¬ 
mize attack effects, and (2) to preserve life and property by operations 
and instructions necessary to reduce attack effects. Civil defense is 
also responsible for supporting services necessary to achieve and maintain 
a capability for effective and coordinated attainment of civil defense 
objectives. It is in this latter category that survival estimating falls. 

Civil defense is the joint responsibility of federal, state, and local 
governments, with the federal government responsible for providing neces¬ 
sary coordination, guidance, and assistance. Federal government agencies, 
including the Office of Civil Defense and the military, advise, guide, and 
assist the states and their subdivisions. Figure 1 shows the major govern¬ 
ment elements likely to be involved in postattack survival estimating. 

The Office of Civil Defense, under control of the Secretary of tic 
Army, is responsible for carrying out DOD civil defense responsibilities, 

OCD has headquarters in the Pentagon and has eight regional offices covering 
the United States and its territorial possessions. The regional offices 
work closely with the state civil defense agencies, tuid through them, with 
local Civil defense organizations. The Office of Civil Defense coordinates 


* "Federal Civil Defense Guide," office of Civil Defense, Washington, D. C. 
March, 1965, 

^ The National Plan for Emergency Preparedness," Office of Emergency 
Planning, Washington, D. C., December 1961. 
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Figure 1 


MAJOR GOVERNMENT ELEMENTS INVOLVED 
IN POSTATTACK SURVIVAL ESTIMATING 



NVCSSC - National Military Command Center System 
Support Center 



























federal activities in support of situation analysis for civil defense 
purposes. OCD also acts us a depository for preattack and postattack 
resource data related to civil defense, and as an exchange point and 
situation analysis center for civil defense planning and operations. The 
National Civil Defense Computer Facility, 0 )>erating under control of 
the Office of Civil Delense, would play an important role in postattack 
survival estimating. 

The Office of Emergency Planning (OEP) is closely associated with the 
Office of the President, and is responsible for determining policy on 
civil defense, and for planning, directing, and coordinating the total 
civil defense program. OEP also determines the civil defense roles of 
other federal agencies and coordinates their civil defense activities with 
each other and with those ol the states. OEP further aids in arranging 
mutual aid compacts among the states and in enacting legislation for civil 
defense purposes. The National Resource Analysis Center (NRAC), a 
facility operated by OEP, coordinates federal activities in support of 
situation analysis for purposes of continuity of government and for central 
programming of resources. NRAC, using its own and associated facilities, 
acts as a depository for preattack and postattack resource data and as an 
information exchange mechanism to support emergency decision making. 

Other federal agencies participate in the national resource evaluation 
programs and provide OCD with data required for civil defense plans, pro¬ 
grams, and operations. They also maintain the capability to assess the 
effects of attack on resources under their cognizance, and provide OCD with 
data required for plans, programs, and operations in support of situation 
analysis for civil defense purposes. Civil defense responsibilities of the 
other federal departments and agencies are presented in some detail in the 
Federal Civil Defense Guide. In general, each federal agency develops 
civil defense plans for use of its personnel, materials, and services 
during a civil defense emergency. Field offices of the federal agencies 
work closely with the states. 

In the event of a national emergency resulting from a nuclear attack 
on the United States, the civil defense mission of the military depart¬ 
ments is to assist civil authorities in restoring order and civil control. 
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by employing whatever military resources are not required for military 
operations. Available military support would be used to return essential 
facilities to operation, prevent unnecessary loss of life, alleviate suf¬ 
fering, and take other actions ns directed to insure national survival 
and a capability on the part of the nation to continue the conflict, Mili¬ 
tary support would be in coordination with, and supplementary to, the capa¬ 
bilities of state and local governments and would be terminated as soon as 
possible, to conserve military resources and to avoid infringement on the 
responsibility and authority of civil government agencies. However, several 
military systems are available for civil defense purposes routinely--for 
example, logistic support by the Defense Supply Agency. The U.S. Army Corps 
of Engineers and the Naval Facilities Engineering Command provide engineer¬ 
ing support in srch civil defense activities as the shelter survey and 
marking program and community shelter planning, as required. A basic 
function of the National Military Command System Support Center, operated 
by the Defense Communicat ions Agency, is damage estimation. This ianlity 
is an important part of the national damage estimation capability, and 
would undoubtedly contribute to postattack survival estimating for civil 
defense purposes. 

State and local governments are responsible for civil defense opera¬ 
tions within their jurisdictions These governments develop the necessary- 
plans, capabilities, and procedures to carry out civil defense activities 
in accordance with state law and federal policies and guidance. State and 
local governments collect data and prepare materials required for estimating 
damage and making situation repor's. Local governments would report infor¬ 
mation to the state, and it in turn would report to the federal agency field 
offices and the OCD regional offices. 
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IV POjjSl BLK DKCISJONS HKD IN TOSTATTACK OPKRAT1ONS 


I! 


It would be futile to try to predict exact postattac'- situations in 
advance, There is no doubt, however, that many important decisions would 
he required in the very early time period following a major nuclear attack 
on the United States, Tills section will identify some of the tyjies of de¬ 
cisions likely to he required and the ]wople who may muke them, This dis¬ 
cussion is centered on decisions to he made at ihe national level. Deci¬ 
sions likely to be required at the local leu-1 will be similar in many 
instances, but their sco|>c and nature may be quite different. 

Some Likely Decisions 

Among some of the more in|>ortant decisions that are likely to be re¬ 
quired at the national level and that will be strongly affected by damage 
assessment results arc decisions regarding: 

• Direction and coordination oi slate and local governments 

• Provision ol area warning, and public information 

• Conservation and distribution oi manpower 

• Relocation ol population to reduce direct threats lo life 

• Allocation of military support to civil authorities 

• Maintenance or reestablishment ol government control 

• Formulation of plans and procedures lor future actions 

• Protection and restoration of essential utilities 

• Allocation of resources and facilities 

• Protection and relocation ol essential stucks and inventories 

« Reallocation or restoration ol communications channels 

• Maintenance or reestablishment of transportation links 

• Relocation of survival resources, such as emergency hospitals 
and medical supplies 

• bstabl1shmenl ol requirements lor additional information 


.fi¬ 
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Early decisions at the lower levels will often be concerned with di¬ 
rect threats to lilo, such as fire and fallout. They will also bo needed 
regarding protection and allocation of food, water, medical supplies, and 
other Inventories and facilities necessary lor protecting life in the 
immediate future. 

The Decision Makers 

The individuals or organisations who make the important decisions at 
the national level may well d«j|iend on who t.he survivors are. Since civil 
defense and national government are an extensive and highly complex net¬ 
work, the decisions to be made will necessarily Involve a large number of 
f Individuals and organizations. Certainly thi' President and his staff will 

| be closely connected with most o! the important decisions, and the Prosi- 

j. dent may well make many ol them hlmsell. The heads oi government depart¬ 

ments will clearly be called on lor many decisions. The head ol the De¬ 
partment of Health, Education and Welfare wou'd unquestionably be called 
i on for decisions regarding the utilization and allocation of medical man- 

4 power and resources. The Secretary of Defense and heads of the military 

| 

departments would tie exacted to assist the President in decisions iegard- 
i ing continuation of the war. The Otiice ol Civil Defense and Office ol 

Emergency Planning would be called on lor decisions regarding: direction 
and coordination ol state and local governments: provision of area wurn- 
| lng and public inlormation; emergency assistance; and many other matters. 

I Although it will often be difficult to identify the decision maker 

in advance, nevertheless th d. Ions must and will be made, and the 
i best possible damage assess"* 1 . formation is needed lor such decisions. 














V POSTATTACK SURVIVAL ESTIVATING PROBLEMS 


The nnturc of t he inputs available) in the* i>om tot tack environment would 
have a significant effect on the character Of the survival estimating sys¬ 
tem. The postal tack actuality would In some ways lie more helpful than pre¬ 
attack assumptions, as in the case of dose rates, where direct reports would 
replace predictions based on assumed yield, average winds, and other fac¬ 
tor a. On the other hand, postattack data could often lie quite misleading 
becuu.se of incompleteness and biases introduced by sensors and communica¬ 
tions. In general, the data following attack are likely to be much more 
varied in content, reliability, and timing than inputs normally assumed in 
preattack planning. 

The survival estimating system must be able to accept these varied, 
biased. and time-phased rejiurts, translate them into a reasonably unbiased 
picture of actual damage conditions, and present findings In a form suit¬ 
able for making operational decisions. Since the data will be basically 
inadequate, it will be necessary to make a variety of assumptions about 
the actual environmental conditions. Systematic rules must be applied 
for acceptance or rejection of these assumptions. Rules must also permit 
modification and repeated testing of assumptions against incoming data. 

Available data augmented by accepted assumptions would provide the rational 
basis for decision. 

[)i (Ter ences between Preattack and Postattack Survival Estimating 

Major differences between the damage assessment problem as it is cur¬ 
rently approached in the preattack period and the problem as it would be 
in an operational, postattack environment have been discussed by Spence 
and Moll.* Some of these differences have been expressed in terms of 
system input and output requirements, and are presented in Table 2. 


* Spence, Richard H. and Kendall D, Moll, "Rapid Damage Assessment in the 
Postattach Environment.” Stanford Research Institute, Menlo Park, Cnltf.. 
February 1961. 
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INPUT-OUTPUT DIFFERENCES iSFTWEEN PREATTACK 
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Source: Spence, Richard H. and Kendall 1). Moll. Rapid Damage Assessment in the Postattack Environment, 

.Stanford Research Institute, Menlo Park, California February 1961. 












Input differences appear to be the most significant. Inputs can gen¬ 
erally be specified exactly in the preattack system, but postattack inputs 
covering actual weapon strikes may t>e approximate, incomplete, or inaccu¬ 
rate. For example, postattack inputs may consist simply of reports of 
nuclear detonations in a general area, whereas preattack inputs provide 
the exact yield and location of each assumed detonation. Wind data pre¬ 
sent similar problems; such data niav be exactly specified for the pre - 
attack case, but there will probably be inaccurate information and uncer¬ 
tainties associated with postattack wind data. For example, wind data are 
provided In general only at G hour intervals, and many changes in wind con¬ 
ditions are possible during such periods. Many attacks no longer assume 
simultaneous bursts, an imp? ivement over the situation that existed when 
Table 2 was pie pared. 

Tiie format lor inputs to the preattack system can be standardized and 
fixed for individual computations. If this were done for the postattack 
system, however, it night well result in unacceptable consumptions of time. 
Thus the postal tack system must be much more llexible than the preattack 
system. Closely related to this h the matter of burst times. Preattack 
systems in the past have often assumed simultaneous bursts, but in post¬ 
attack systems, burst times will probably occur and t>e reported at varying 
intervals. These and other sporadic inputs must be processed as they be¬ 
come available, and then revised as new information is received. Stag¬ 
gered reports will be received from vaiving sources, in contrast lu pie- 
attack analyses assuming a single group ol dala from one source. 

Preattack systems are generally designed to produce the total output 
in a single computer run from one set ol input conditions. In the post¬ 
al tsek period, on the other hand, decisions to be made require that output 
be produced on request by users, or whenever damage information changes. 
Postal tack decision makers will probably want periodic reports as well as 
the ability to interrogate the system on specilic i*oints oi interest. 

Formats for preattack systems are generally designed to minimize pro¬ 
gramming and computer computation rather than to provide ease in interpre¬ 
tation and analysis oi output. For a postattack system, easily read for¬ 
mats based on rapid analysis in an operational environment are necessary. 
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VI OBJECTIVES OF TTvANSATTACK AND POSTATTACK 
SURVIVAL ESTIMATING 


The traditional role of survival estimating systems is to estimate 
the losses of population and resources resulting from nuclear attacks. 

As survival estimating techniques have become better understood and de¬ 
veloped, it has become apparent that the basic techniques can also be used 
to synthesize an attack during the transattack and early postattack pe¬ 
riods. This latter capability could be very useful as the basis for warn¬ 
ing populations and taking countermeasures. 

Operational survival estimating systems are similar in many ways to 
preattack planning systems, but significant differences exist. For exam¬ 
ple, an operational system of survival estimating should provide for veri¬ 
fication and correction of indirect damage assessments. The range and 
variety of attack parameters considered in an operational system will prob¬ 
ably be more limited than in preattack planning systems. Furthermore, the 
resource base upon which the damage assessment is made should be formulated 
from larger groupings, i.e. , the unit size of the data base can be larger, 
or the data base can be sampled. 

The attack synthesis function of operational survival estimating per¬ 
forms two important functions: (1) it permits the extension of resource 
damage analyses into nonreporting areas, and (2) it can be used to estimate 
the strategy being employed by the attacker.* It is expected that the 
weapon data needed for damage and resource evaluation will be available 
from automatic and manual sensing systems, and from direct reports based on 
aerial and ground observations. Since portions of the sensing and intel¬ 
ligence gathering network are likely to be destroyed during the attacks, the 
nonavailability of complete information must somehow be compensated for if 


* This in turn can be used to predict the development of the attack. 


* 










the damage estimates are to be complete. One way to accomplish this pur¬ 
pose is to employ a combination of judgment and inference from available 
information relying mainiy on the latter. Synthesis of the attack strategy 
could also be useful in the early stages in estimating whether the attack 
is oriented toward civil population, military retaliatory forces, indus¬ 
trial resources, or other objectives. This type of attack-pattern projec¬ 
tion and its ability to identify probable targets as the attpck progresses 
could be crucial to population or other resource warning systems. 
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VII A STATISTICAL CONCEPT OF SURVIVAL ESTIMATING 


| 


General 

The requirements to be placed on an cqierat j ona 1 survival estimating 
system, together with the probability of limited available information, 
suggest a statistical rather than deterministic approach. Concepts and 
some details for a statistical approach are presented in this section. 

The presentation of this section is largely in terms of generalized flow 
diagrams with accompanying descriptive material. More detailed flow dia¬ 
grams are given in Appendix A. 

The system described is Intended to operate with incomplete and per¬ 
haps inaccurate data. The underlying principle is that of statistical 
inference, i.e., the definition of an entire attack design from an early 
and probably small sample of attack reports. The statistical inference 
techniques devised are then coupled with heuristic reasoning to provide 
the information required to implement the survival estimation procedure. 

Unique features associated with the survival estimating sjstem pre¬ 
sented include: 

1. The system operates with incomplete and inaccurate data, i.e., 
only a sample ol the data describing the attack is needed to 
produce preliminary survival estimates. The approach in this 
respect is therefore somewhat similar to that followed by the 
major radio and television networks on election nights in which 
the networks attempt to estimate the final vote from preliminary, 
and other scanty, information. 

2. The system makes extensive use of direct damage reports. Di¬ 
rect damage reports are used to infer NUDET-type information 
and to verify other reports. 

3. The system identifies requirements for status reports that are 
needed to upgrade system results; that is, the system is designed 
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to automatically generate requests for data where there are gaps 
In the data, or where necessary data have not been provided. 

4. The conceptual system described by this report is a real time sy^ 
tem in that input data are continuously processed and handled 
statistically to infer attack size and target system. Production 
of survival estimates, the next step, can be essentially real 
time if the parametric survival estimating method is used. How¬ 
ever, generation of survival estimates using large data bases 
would require appreciable time periods, thus falling outside 

the usual definition of real time. Also, it should be realized 
that detonation and direct damage reports will in most cases not 
be received as soon as detonations occur; there will be some de¬ 
lay before this input data enters the system. 

5. The system may also identify certain target areas in advance of 
their being attacked, based upon the attack pattern produced from 
the sample of attack data. 

6. The synthesized target design can be coupled with a variety of 
survival estimating techniques depending upon the needs and de¬ 
sires of the system operator. 

Sta tistical Damage Assessment 

A highly relevant statistical method of estimating damage from nuclear 
attacks was devised and presented by Laurlno.* The general logic of the 
overall damage assessment procedure presented by Laurlno has been summar¬ 
ized by Bothun and is shown in Figure 2. The system logic as presented 
by Laurlno consisted of seven distinct phases. 

1. Screening of incoming reports 

2. Definition of the target system 

3. Determination of average yields per target 

4. Determination of fraction surface yield 

* Laurlno, Richard, David Goodrich, and Donald Doane, "Statistical Meth¬ 
ods of Estimating Damage from Nuclear Attacks on CONUS (U)," Stanford 
Research Institute, Menlo Park, California, 1962. (SECRET) 
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Figure 2 

SUMMARIZED LOGIC OF DAMAGE ASSESSMENT PROCEDURES 



Source: Bothun, Ft. B. , Uses of Radiological Fallout Information in Operational 

Pos tat tack Damage Assessment Systems (U)," Stanford Research Institute, Menlo Park, 
California, January *965 (CONFIDENTIAL) 
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5. Augmentation of target system 

6. Determination of attack parameters 

7. Determination of damage estimates 

The proc',dures presented by Laurino provided gross estimates of popu¬ 
lation casualties and resource damage within the limitations of the data 
base and the approximating procedures used. The procedures provided only 
for survival estimates derived from precalculated survival functions of 
yield versus damage. In the event of computer malfunction, and where 
precise location of damage is not required, this system makes possible 
quite rapid manual calculations that take into account the pattern and 
weight of attack. This feature is carried over as a capability of the 
system presented by this study. Laurino's system is based on a minimum 
number of independent parameters, and the submodels of the system are as 
simple as possible, consistent with the general requirements for accuracy 
and flexibility. 

Many of the statistical concepts presented by Laurino are employed in 
the overall system devised in this study. Where Laurino's original con¬ 
cepts have been adopted, they have been further developed and extended. 
However, several major differences exist between this and the Laurino 
study, including: (1) provision is made lor extensive use of direct dam¬ 
age reports; (2) the system identifies and generates requirements for 
status reports; (3) the system is designed to ojterate on a real time basis, 
(4) techniques are presented coupling the basic target identification sys¬ 
tem with several survival estimating techniques rather than one; and (5) the 
statistical methodology and techniques have been modified and improved. 

Attack Synthesis Coupling with Survival Estimating Systems 

A brief description of various survival estimating and damage assess¬ 
ment systems will be presented before proceeding to the specific subsystem 
considered in this study. The purpose of this discussion is to summarize 
the various survival estimating and damage assessment systems that might 
be employed with the statistical target synthesis concept of this report. 

The statistical target design system produces NL'DET-type information 
for each actual ground zero determined, This information includes weapon 
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yield, Height ol burst, ground zero, und time of burst. The data rnav be 
summarized for use with parametric systems (based on preattack planning 
calculations) or may tie used with individual ground zeros for more accu¬ 
rate survival estimating systems. 

Some of the survival estimating and damage estimating systems that 
might be employed are shown as a function of time in Figure 3, The time 
axis of this figure is not intended to lie exact, but only somewhat repre¬ 
sentative. The time (leriml shown ranges from tl-lmur (considered to be the 
time of the first detonation), to the time when reporting is completed, or 
essentially completed, 1rora all attested areas. 

The first system indicated in Figure 3 is that which employs preat¬ 
tack calculated damage functions. The principal advantage ol the para¬ 
metric system is that it is very rapid and car, easily he accomplished manu¬ 
ally if necessary. The principal disadvantage ol this system lies in the 
uncertainty regarding its accuracy. Unfortunately, the attack designs se¬ 
lected tor the precalculated damage (unctions may not closely resemble the 
actual attack. The achievement of surprise is a fundamental military 
principle, and history is replete with examples in which surprise has been 
achieved. Thus, an actual attack on this country may not compare well with 
assumptions made for any of the preattack calculated damage lunctions. 

The second major category of sur\i\al estimating and damage assessment 
systems that might be employed is based on a combination of weapons ef¬ 
fects scaling models and resource data bases. Two such systems are illus¬ 
trated in Figure 3--tho.se based on .SKI SAMPKO ami DASTAP programs. These 
systems are currently used for preattack planning studies, but they may 
easily be adapted lur survival estimating in transattack and postal tack 
jiericxis. liASTAP employs a rather complete data base ol some 44,000 re¬ 
source points. The SAMPKO program employs only a small fraction of 
IMSTAP's resource points—about 4,000. Thus the SAMPKO program may be 
run with much less computer time than the complete data base DASTAP pro¬ 
gram. The running of a SAMPKO program may require about 5 minutes or 
less, whereas the running ol a DASTAP program may require an hour or more. 
Although SAMPKO runs are less accurate than DASTAP runs, SAMPKO is accu¬ 
rate enougti to be suitable for quick runs concerning national estimates. 
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.SAM PRO thus appears to lx- eminently suited for use with statistical 
systems, when the sample is relatively small. The principal disadvantage 
of the SAMPRO system as compared with the parametrie system is that SAMPRO 
requires a computer, and the parametric system may be o|>eratod manually. 
However, this advantage may be more theoretical than actual, because incoming 
reports from affected areas are likely to be so numerous, jurticularly 
those concerning radiological ha/.ard conditions, that a computer will be 
needed for processing the quantity of incoming data in any case. 

The next survival estimating system depicted In 1’lgure 3 is one based 
on a combination oi weapons effects scaling and reports 1 rum the field. 

This system would employ the statistical target synthesis method presented 
by this report, together with weapons eilects staling models m the very- 
early transattack and pustatiack periods. In this respect, it would be 
like the SAMPRO and DASTAP based survival estimating models described 
above, but would also diller I roi,. them In using direct reports I rom af- 
lected areas to correct and update quick survival estimates made in the 
very early periods with .SAMPRO or DA.3TAP. That is direct damage and 
status report.- iront the various resource |joints in the livid will be used 
to provide new survival probabilities at each re ouree point and replace 
earlier survival probabilities derived 1rom weapons effects scaling models. 
Thus, in this particular model, survival probabilities at a given resource 
point will initially be based only upon the statistical target identifica¬ 
tion system and weapons oflects scaling. These survival probabilities will 
all be replaced as direct damage reports and status reports are received 
from the oiierating areas. A transition period would exist, in which the 
total survival est: ate would be partly based on weapons eilects scaling 
and partly on direct reports. This system has the advantage that it can 
provide a continuous total nit tonal survival estimate using the best 
available data at all times. This model has not been developed, anil 
considerable additional research will be required to complete it. 

The remaining survival estimating system indicated in figure 3 is one 
based on direct reports only. '’Tie obvious disadvantage oi this system is 
that it provides no national survival estimate in the transattack and 
early postattack jx*nods. Smli a system is clearly not suited lor meeting 
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civil defense objectives. Furthermore, it has no advantage over the sys¬ 
tem described in the paragraph above. 

The preceding discussion has described how a statistical target syn¬ 
thesis system might be used with various survival estimating and damage 
assessment systems. In general, it appears that the capability for carrying 
out a combination of these systems is desirable. The parametric system may 
be useful for making flash estimates in the very early periods. If the 
parametric systc;.. is used, it should be followed closely by the SAMPRO sys¬ 
tem which can provide improved estimates in very short times. The most 
accurate and useful system in the late - periods will be a system like that 
described above which incorporates the best features of weapons effects 
scaling models and direct report models. All of these systems should be 
useful at the national level. 


Target Information 


Certain target categorizations are desirable in the development of 
the statistical concept. First of all, it is assumed that the operator 
of the system has available to him a reasonably complete and well-defined 
target fist. This target list will presumably be similar, or identical, 
to target lists currently used in preattack planning studies. 

The target list will be divided into target categories, and the total 
target list might typi«' ally cor. In some 20 to 30 target categories. For 
example, a target category might consist of Mlnuteman sites. Other ex¬ 
amples of possible target categories include submarine bases, Strategic 
Air Command bases, and population centers. 

For purposes of developing decision rules and estimating uncertain 
weapon parameters, the target list, will also be broken down by vulnera¬ 
bility groups. The three vulnerability groups used are designated as 
military hard, military soft, and nonmilitory. The ways lr, which target 
categories and vu lerabllity groups are used will be described in follow¬ 
ing sections. 


The statistical large! synthesis system also uses a 
list, by category, provided that such a list is 


target priority 

The availability 


avallable . 












of a target priority list is helpful but not essential for the system. 

Such a list would presumably be based in large part on strategic intelli¬ 
gence 2 nd pattern analysis of the early attack reports. It would also 
presumably be based upon the assumption of a rational enemy, but not nec¬ 
essarily. Use ot a target priority list will be further explained later 
in this report. 

Components of the Statistical Target Symnesis System 

The statistical target synthesis system is made up from a number of 
subsystems. The subsystems include data input and verification, direct 
damage report processing, requested status report processing, detonation 
record maintenance, initial target selection, and final target system 
definition. Each ol these subsystems will be described briefly below. 

Tin inputs, principal steps to be accomplished and outputs of each sub¬ 
system are shown in associated figures. A greater level of detail for 
some subsystems is indicated in Appendix A. 

Data Input and Verification 

The principal elements of the data input and verilication subsystem 
are summari ■ d in Figure 4. 

Subsystem Inputs 

The subsystem is designed to receive inputs irom many different 
sources. The principal sources are thought to be sensor NUDET, manual 
NUDET, satellite, direct damage reports, and requested status reports. 

The receipt and initial processing of this information presents many var¬ 
ied problems, Some of the data received will probably be accurate and 
some of the data will, no doubt, be highly inaccurate. The report reli¬ 
ability from the several sources may vary greatly. In the extremely early 
periods (e.g., the first detonation), there will no doubt be insufficient 
data for making reliable estimates. As time progresses and even in reas¬ 
onably early periods, vast quantities of data may become available. There 
are many thousands of radiological teporling stations, and many of these 
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stations will be reporting early. Dozens of manual NUDET reports cov¬ 
ering a single detonation may be received. In this kind of situation, 
a system not capable of handling large quantities of data, and selecting 
important data while rejecting less important or trivial data, would be 
completely ineffective, and possibly inoperative. 

The subsystem assumes the availability of electronic sensors for re¬ 
ceiving and reporting essential NUDET data. If such a system is actually 
available during an attack, it presumably will provide weapon yield, height 
of burst, ground zero, and time of burst. The accuracy of each piece of 
Information may vary considerably, depending on such things as distance 
between sensor and detonation, and other factors. Some elements of the 
data will probably be more accurate than others. One might expect the 
time of burst to be quite accurate, whereas the reported ground zero might 
not be very accurate if the sensor and detonation are separated by consid¬ 
erable distance. 

Manual NUDET reports made by trained and untrained observers will 
probably also be received. Various techniques for estimating NUDET infor¬ 
mation manually are given in apprcpriate civil defense and other govern¬ 
mental publications. For example, weapon yield may be estimated from 
cloud height. The distance separating the observer and ground zero may 
be estimated by the time interval between the initial weapon ilash and 
the arrival of the blast wave at the observer's location. These and other 
techniques may be used to estimate NUDEr data manually. In general, one 
would expect the accuracy of manual NUDET data to be considerably less than 
that of Sensor NUDET data. However, the existence of one or more manual 
NUDET reports for a given detonation would enable one to accept with some 
confidence that a detonation had actually occurred. Faulty or incorrect 
data concerning the existence of a detonation might conceivably be re¬ 
ceived from electronic sensor systems. 
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A further main source of data indicated in Figure 4 is direct damage 
reports. These reports will be received in accordance with various civil 
defense directives.* Such reports may be received in the form of NUDET 
data, radiological contamination data, or other ways, such as a simple 
description of damage sustained. Direct damage reports will also pre¬ 
sumably be received from civil defense operating areas, as in the manner 
described in certain current civil defense documents.^ In this concept, 
proposed by Strope, each civil defense operating area would report its 
status by specified basic operating situations. The categories proposed 
give basic data regarding the fire and fallout condition of the operating 
area. This concept will be discussed further in connection with the 
direct damage report processing subsystem. 

As Indicated previously, the statistical target synthesis system 
will generate its own requests for operating area status when initiai 
data are needed. Status reports would generally take the form of direct 
damage estimates and would be converted to a similar format when received. 

The subsystem is designed to operate in the absence of one or sev¬ 
eral of the above sources of data. Presumably, direct damage reports and 
requested status reports could be made available In any case, but perhaps 
not for all operating areas. An essential point is that the subsystem is 


* ".National Warning System (NAWAS) Operations Mr nua 1," Department of 
Defense, Office of Civil Defense, iJecember J9J6. 
f "Concept of Operations Under .Nuclear Attack,' 01 I ice of Civil Defense, 
Washington. D.C., June 26, 1967. (Working Draft) 
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flexible enough to accept inputs from many different sources and is not 
dependent on any single source lor its operation. 

Subsystem Accomplishments 

Data coming from the attack environment will first be tested for con¬ 
sistency. This test is intended to exclude data which are clearly inaccu¬ 
rate. For example, a consistency test might be one in which it is ascer¬ 
tained whether the reported ground zero is within range of the sensor. 
Another possible consistency test might be to ask whether the reported 
time of burst is later than the last O.K. (unaffected or undamaged) status 
report. Still another possible consistency test might have to do with 
weapon yield. Strategic intelligence should tell us what range of yields 
in general might be anticipated. If a weapon yield outside these limits 
is reported, it should be subjected to scrutiny by the system operator, 
and if consistency is questionable, status reports from these particular 
ojierating areas should be requested (as for other questionable cases). 

After consistency tests, the next activity will be to transcribe 
incoming reports to punched card form in a standard format. Once this has 
been done, these reports are referred m as detonation records. Detona¬ 
tion records aiv initially derived only Iron sensor NTDET, manual NUDET, 
and satellite data. The elements ol data that are proposed for each det¬ 
onation record are as iollows: 

1. Location and type of rejxirTing unit 

2. Yield of weapon 

3. Height of burst (HG3) 

4. Time of burst (TOD) 

5. Ground zero (UTM coordinates) 

6. Report number lor later identification 

Subsystem Outputs 

The outputs irom the data input and verification subsystem are det¬ 
onation records, direct damage records, and requested status records. 
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These outputs will presumably be in the form of punched cards. They will 
be immediately available as inputs for other subsystems. 

Possible Future Extensions of the Subsystem 

The subsystem activities described above are by no means intended to 
be a complete listing. As new sources of data become available, they 
should be considered for incorporation into the subsystem, For example, 
it has been reported that satellites will provide photographic data. These 
data should be used to obtain or supplement detonation records. Various 
Information sources may also provide physical census data which can be sum¬ 
marized or reduced to obtain detailed specific results of an attack.. 

Direct D amage Report Processing Subsystem 

The principal inputs, accomplishments, and outputs of this subsystem 
are indicated in Figure 5. 

Direct damage records are treated somewhat differently in the data 
input, and verification phase since they may vary considerably and probably 
will not contain all of the elements of data required to produce a com¬ 
plete detonation record. Direct damage records should provide, as a mini¬ 
mum, the location of the operating area, identification code of operating 
area, the status code, and time of status. Various status codes might be 
employed but the status code considered here is that presented in "Con¬ 
cept of Operations Under Nuclear Attack," Office of Civil Defense, Wash¬ 
ington, D.C. The basic operating situations and their corresponding codes 
are indicated in Appendix B, For example, a heavily damaged area might 
carry the following code designation: 

Basic 

Map Operating 

Code Situation Situation Definition 

b L0RAD-HIF1RE Dose rate between 0.5 and 50 r/hr. Many 

fires beyond control capability. 
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SMSAs or other geographical areas. 

Does this damage record satisfy a request for status? 















This particular code will probably change as new reporting systems 
develop, but it is representative of possible simplified damage codes. 
Requested status records should indicate facility and location code, 
status code, and time of report. 

System Inputs 


The inputs for this subsystem are the direct damage file, direct dam¬ 
age records, and requested status file, Direct damage records are received 
from the data input and verification subsystem and are sorted by location. 
Direct damage records also receive a minor or secondary sort by time of 
report. The direct damage file is created upon receipt of the first di¬ 
rect damage record, and is maintained thereafter. 

Subsystem Accomplishments 

When a new direct damage record is received, the iirst action is to 
update the direct damage file. This is done by changing the status of op¬ 
erating areas and facilities that have previous reports so that the status 
agrees with the new report. If no rejjort had previously been received from 
the relevant operating area or facility, the status would he changed from 
"no report" to the current status code. 

One of the outputs of this system is a change and summary report which 
is presented to the system operator or other interested persons for system 
review. The next activity ol the system, therefore, is to record all 
changes and operating area and facility status for the change and summary 
report. 

A further activity of this subsystem is to provide indicators to the 
system operator regarding the overall status of the reporting system. Vari¬ 
ous indicators might be employed --1 or example, the ratio of the number of 
reporting operating areas to the total number of operating areas in the 
system would be of interest to the system operator. Another possible 
indicator would he the fraction of the system not reporting by areas. 

Such indicators should provide a measure of overall performance oi the 
reporting system. 
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One of the principal uses of direct damage reports Is to provide data 
for generation of detonation records. This is discussed in the section on 
detonation record maintenance and in Appendix B. 

Tile iinai step in direct damage report processing is to determine 
if requested status records can be satisfied by a search of records In 
the direct damage file; that is, it should be ascertained whether the 
incoming direct damage record can satisfy a current status request. 

Subsystem Outputs 


Outputs of the direct damage reporting subsystem include an updated 
direct damage file, U|xiated requested status file, and change and summary 
reports for system review, 

Requested Status Report Processing Subsystem 


Inputs, activities, and outputs of this subsystem are indicated in 
Figure 6. In general, the handling of requested status reports is sim¬ 
ilar to that given for direct damage reports. 

System Inputs 

System inputs for the requested status report processing subsystem 
are the requested status records themselves, the requested status file, 
and the direct damage file. The requested status records should be sorted 
by location with a minor sort by time of report. 

Subsystem Accomplishments 


Status reports will result only in response to a s(>ecific request. 
Therefore, one of the first and [perhaps most important actions will be to 
update direct damage file records for appropriate operating areas and 
fact 1ities. 

The requesting system should also tie cleared of requests for status 
that have been received. Finally, reports o! changes in detonation rec¬ 
ords ami direct damage files should be printed out lor inspection bv the 
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system operator aiul other interested persons. As with direct damage proc¬ 
essing. detonation records util ultimately be inferred 1rom requested 
status records as needed. 

■Subsystem Outputs 

The principal output ol this subsystem Is an updated direct damage 
file. Other outputs are an updated requested status tile, and printouts 
of changes in files and status of lacililies cunfirmed by requested status 
records. 

intonation Record Maintenance Subsystem 

Some of the statistical procedures covered n. this subsystem are 
described in greater detail m the next chapter ol this report. Several 
of the key points relating to the statistical basis lor the survival esti¬ 
mating concept center around tins part cular subsystem. The principal 
components of this subsystem are indicated in Figure 7. 

Subsystem Inputs 

The inputs lor tills suhsvsten are detonation records, the detonation 
record file, the direct damage file, and the estimated actual ground zero 
file. Detonation records arc soiled by nmc ot burst and by geographical 
1 cent i(uis. 


Subsystem Accomplishments 

As new detonation records are received, they are used to update the 
detonation :ccord lilc. Su|>erseded records should be replaced, and new 
records should be added to the tile. In addition, detonation records 
should be compared with relevant direct damage records, to determine 
whether they are compatible. If the tie tonal ion record is not compatible 
with actual status, then the- detonation record should tie modified to 
iliiiunate tile incompatibility. 





detonation record maintenance 



So r t ed by 










The number of bursts represented by the sorted Jetun.it inn records and 
necessary we a port parameters are next determined. Tins • s dune statisti¬ 
cally, as described in the next chapter. The number ol bursts is estimated 
by an algorithm that separates the* information inputs into g • < h; depend¬ 

ing on the overlaps of confidence regions associated with tin fts.ritl 
ground zeros and times of burst. bach group generated tn ssr 4 sisa 
represents one burst. 

Detonation records are grouped by estimated actual ground r.v tot. 

These ground zeros are determined statistically from the azi.mul ti data by 
a least-squares method that is described 111 tne nc-xt chapter ol this re¬ 
port. Standard deviations o! the relevant reporting media are u c c-u in con¬ 
nection with determining location of reported ground zero, weapon yield, 
ami time ol hurst. 

The next ac.ivity ol the subsystem, and a key activity lor operation 
ol the whole system, is the generation ol estimated actual ground zeros 
and related weapon jiaramct er.s 1 rum direct damage records and status re¬ 
ports, The generation ol these data 1 rom (Intel damage reports and status 
reports is explained in greater detail in Apgietulix B. In general, detona¬ 
tion' are Inferred by pattern recognition techniques. A detonation could 
presumably be interred Irom a single direct damage lecord. or status re¬ 
port, but the availability ol mure than one iecotu should greatly enhance 
Lhe accuracy and value ol detonations interred 1 'trough pattern recognition. 
As time passes and none status reports are received, new estimated actual 
ground zeros may be generated lor the same i'eI onelion, Thus, care should 
be taken to update the estimated actual ground zero file. 

Subsystem Outputs 

The principal output ol this subsystem is a description ••• the attack 
size. It yields the total number ot estimated actual ground zeios and re¬ 
lated detonation record information, determined on a statistical basis. 

The target categories that have been attacked may lx? obtained 11 desired. 
Subsystem outputs also include change reports and an updated detonation 
record file-. 
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Initial Target Selection Sub system 

The principal components of the initial target selection subsystem 
are indicated in Figure 8. The main purpose of this subsystem is to make 
initial target selections from the various target candidates. 

Subsystem Inputs 


One of the principal inputs for this subsystem is that derived from 
the detonatton record maintenance subsystem--i.e., a file of estimated 
actual ground zeros together with associated weapon parameters (NUDET- 
t\|>e information) as well as target categories that are associated with 
estimated actual ground zeros. Other injjuts for this subsystem include 
the direct damage file, detonation record file, and the target candidate 
file. 


Subsystem Accomplishments 

The first principal action of this subsystem is to estimate the feasi¬ 
bility of estimated actual ground zeros by checking the status of facili¬ 
ties in the area. The detonation is highly sus|>ect and should be removed 
from the estimated actual ground zero 1i le li either ol two conditions are 
met: if there are facilities with "no hit" status reports latei than the 
estimated time of burst; or 11 the reported damage is grossly incompatible 
with the estimated characteristics ol the detonation (e.g., ail facilities 
in t ne area rc*i>ort light blast damage and no significant fallout, while 
the detonation is supposed to have been a 10 MT surface burst). Also, 
infeasible detonations are primed out lor review by the system operator. 

Each detonation record group (associated with an estimated actual 
ground zero) is then compared against a list of target candidates. The 
hypothesis that a given target candidal-a was the target lor the group is 
tested bv u*-ing rang* and azimuth data to compute a variate that is dis¬ 
tributed approximately ihi-square , When the chi-square test fails to 
yield a possible target (l.e.. a given detonation record group apparently 
has no feasible target), this information is printed out for the system 
ope ru tor. 


do 



























The next principal action is to check all possible targets for con¬ 
sistency with the associated detonation parameters, i.e., yield, height 
of burst, etc. It may turn out that there are no possible targets con¬ 
sistent with the associated detonation (i.e., the estimated weapon charac¬ 
teristics are not appropriate for use against the vulnerability group or 
groups of the possible targets). In this case, this information is printed 
out for the system operator. Heuristic methods are then used to make the 
final selection of a target (or targets), if more than one target candidate 
passes the chi-square and consistency tests. 

For example, a likelihood function might be used to favor the selec¬ 
tion of target candidates in those target categories that have a higher 
percentage of reported hits. The reliability levels of associated reports 
might also be used. 

Subsystem Outputs 

The principal output of this subsystem is an initial target list. 

The various environmental Inputs of the overall system, including NUDET 
reports, dire< damage reports, etc., have been analyzed and processed 
statistically to produce an initial list of targets based on a sample of 
the environment. The subsystem also prints out a listing of estimated ac¬ 
tual ground zeros with no selected targets, and a list of infeasible esti¬ 
mated actual ground zeros. 

Final Target System Dellnltlon Subsystem 

This subsystem takes the initial target list and develops an estimate 
of the complete target system, which in turn becomes the basis for sur¬ 
vival estimates; that is, in this subsystem the full attack design is 
estimated from a sample of the targets that have presumably been attacked. 
The principal elements of this subsystem are indicated in Figure 9. 

Subsystem Inputs 

Subsystem inputs include the initial target list, the direct damage 
file, and the target candidate file. Priorities are a-signed to the 
target candidate- by category. 
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TARGET SYSTEM DEFINITION 
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Subsystem Accomplishments 


The principal activity of this subsystem is to develop the full at¬ 
tack design. This is done by accepting as part of the final target system 
all of the target cagegories that have sustained hits, and by selecting 
additions from the target priority list. All targets in each target cate¬ 
gory that have been hit are included on the final target list unless no-hit 
reports have been made for that particular target candidate, with time of 
report later than the relevant time of burst. 

Target categories are considered for inclusion in the final target 
system ever, though none of the target candidates in that category have 
been reported hit, provided the target category is higher on the priority 
list than accepted target categories with no hits, or is at the same level 
with them. However, decision rules are applied to target categories thus 
added to the list, and such categories are rejected when no other cate¬ 
gories in the associated vulnerability group have been hit or when there 
are no-hit status reports for target candidates in that particular target 
category. Status reports are requested on all unreported targets added 
to the target system if status reports have not already been requested. 
These targets should be warned that they may be on the enemy target list. 

Certain weapon yield determinations must then be made in connection 
with target candidates included in the final target system for which 
weapon yields have not been previously determined. Average weapon yields 
and average fraction surface yields arc determined for each target cate¬ 
gory and for each vulnerability group. Appropriate decision rules are 
applied so that needed weapon yield information for the final target sys¬ 
tem can be determined. See Appendix A. 

When no ground zero is available for a target that has been added to 
the target list by i iference, the target location is used as the esti¬ 
mated ground zero. 
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Subsystem Outputs 


The principal output of this subsystem is all of the information that 
is needed as input information for paramet-ic DASTAP and SAMPRO type sur¬ 
vival estimates. For the parametric system, the essential inputs are total 
and surface yields by target category. NUDFT-type information is required 
for the DASTAP and SAMPRO systems. The subsystem also provides an updated 
requested status file and a list of target candidates to be warned. 


Sample Target System 


Before proceeding with a description of how the final target systv.ni, 
as described by the preceding sections, is coupled to various survival 
estimating systems, an illustrative example is given of how the final tar¬ 
get system selection process operates. A simplified target system, con¬ 
sisting of only fou. target categories and up to six targets j>er category, 
is shown in Figure 10. For purposes of this example, it is considered 
that the target candidates can be placed in one of four classes as follows 

1. On the initial target list ana having confirmed damage 

2. On the Initial target list, but having no confirmed damage as yet 

3. Confirmed "no damage" 

4. No report at all. 

Other target codes, such as damaged but not selected for the initial tar¬ 
get list," are possible and might be included. 

Target category 1 contains an example ul each of these classes. Tar¬ 
gets 1 and 3 are or. the initial target li^t, and damage to these targets 


has been confirmed. 


vrgel 2 of this category is on the initial target 


list, but damage has not yet been confirmed. No report has been received 
from target 4, and a confirmed "no damage" report has been received from 
target 5. Target category 1 is thus selected for inclusion in the final 
target system, except that target 5 is deleted. Status reports will have 
been requested for targets 2 and 4, where no confirmed damage ana no re¬ 
ports have been received, respectively. 
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SAMPLE TARGET SYSTEM 
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Target categories 2 and 3 provide further examples similar to those 
given for target category 1. Target category 4 is not selected for inclu¬ 
sion in tlie iinal target system since no targets in that category have been 
placed on the initial target list. .Similarly, status reports are not re¬ 
quested for the target candidates ui target category 4 lor which there are 
no reports, since there is no apparent activity in this target category. 

If target categoty 5 is equal to or higher than target categories 1, 2, 
or 3, it would also be 1nc Luded in the final target system. This would 
depend on its meeting other decision rule criteria that may be established. 

The simplified example described above applies to the final target sys¬ 
tem at some specilied point in time. In actual practice, the definition of 
the final target system, being a dynamic process, is constantly improving 
with time. 

Parametric Survival Estimating 

A parametric subsystem lm sun iv.il estimating, with l n 1 orma t i on pro¬ 
vided by the Iinal target system del inition subsystem, is indicated in Fig¬ 
ure 11. Tins is a high is simplified s s stem that could be either manually 
or machine ojie rated. 

Subsys tem Inputs 

Inputs for this subsystem include total and surface yield ior the at¬ 
tach by target category and precalculated survival functions. The survival 
1 unctions will be based on preattack planning calculations for various 
types and weights of attack. Such survival functions are currently avail¬ 
able for a wide range of ty|jcs ; of attack. One example of this kind of sur¬ 
vival Junction is given, in Figure 12, "Expected Fatalities Versus Weapons 
Detonating at U.S. Cities."* Figure 12 enables one to make gross predic¬ 
tions of the percentage survivors as a j unction '>! total attack design. 
Important parameters to be considered in such functions include attack ob¬ 
jective, height of burst, and shelter posture. The assumption- made will 
have important eflects on the lunctions. Similar survival functions are 
available for survival estimates concerning damage to various types of in¬ 
dustry ami other resoui ces . 


* Exemplar only. m>t based on actual calculations. 
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As suggested earlier in this report, however, there is no assurance 
that the actual attack will closely resemble the hypothetical attacks upon 
which precalculated survival functions are based. This is a major limi¬ 
tation of parametric survival estimating based on precalculated survival 
functions. Such functions art c. urly applicable only to na nal esti¬ 
mates. For example, the pre ice or absence of ballistic r. ssne defense 
systems could have ir,. ortant i plications for precalculated functions, as 
well as for the planning parameters mentioned above. 

Subsystem Accomplishment and Out put 


The principal function of this subsystem is to provide quick national 
estimates of surviving population and other resources. Steps required to 
carry out this function are indicated in Figure 1], 

DASTAP and SAMPRO Survival Estimating Subsystem 

Survival estimating based upon the DASTAP and SAMPP.O models is suit¬ 
able only for machine operation. The DASTAP model calculates survival 
probabilities, standard intensity, and equivalent residual dose for each 
of some 40,000 resource points (depending on the data base used). Results 
of these calculations are inputs to an environment tape which, in turn, 
is used to develop total survival estimates. Other inputs to the environ¬ 
ment tape include wind patterns, weapon descriptions, damage iunctions, 
mean lethal radii, and yield scaling features. Coiru>arable survival models 
are in use by NRAC, RAND, and other organizations. The principal limita¬ 
tion of this type of model is the time required for running the program. 
Running times range typically from 1 to 5 hours, depending upon the type 
ol attack and whether or not fallout is present. Attacks involving fall¬ 
out require more time. 

The SAMPRO model is a modification of DASTAP and employs statistical 
sampling techniques to reduce the number of required resource points to 
about one-tenth ol those used by the DASTAP model. Running time is thus 
greatly reduced and there is very little sacrifice in accuracy for na¬ 
tional estimates. In general, the SAMPRO model appears to he quite 
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suitable for use witlx a statistically designed target estimating system. 
The principal component s of the DASTAP and S.Uil’Ki) survival estimating si' 1 '- 
systc„.s are indicated in Figure 13. 

Subsystem Inputs 


Inputs for the DASTAP and SA.MPHQ models include: (1) resource point 

aata, (2) weapon data, (3) wind data, (l) damage functions, and (5) shel¬ 
ter conditions. Weapon data will be derived from the final target system 
definition subsystem. Wind data and shelter conditions may lie assumed cr 
otherwise obtained as desired. 

Subsystem Accomplishment and Out put 

The DASTAP and SA.MPHO models produce environment tapes ol survival 
probabilities, H + 1 intensities, equivalent residual dose, and other 
parameters. The enviroiiment tape is then edited against resource points. 
This edit results in survival estimates for population and other resources. 

Survival Estimating Systems 

This sect ion has presented methodology (or estimating the magnitude 
and scope ol nucleai attacks, and lor combining this methodology with 
existing damage assessment techniques li> produce national survival esti¬ 
mates. Thus the objectives of the study are mot and exceeded, since more 
than one survival estimating system, each capable of coping with, incom¬ 
plete attack information, has been devised. In the opinion ol the auth¬ 
ors, systems based on the parametric and SAMPItO methodologies should be 
fully implemented lor use by the National Civil lK.fer.se Computer Facility. 
Further, another system that should be fully developed and implemented is 
a survival estimating system based on the statistical methodology lor es¬ 
timating attack size and damage assessment techniques that collate weapons 
ellects with lacilily status and damage reports. The next section gives 
details regarding, statistical methodology lor estimating attack size. 
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VIII STATISTICAL METHODOLOGY* 

Genera 1 

This .section of the > o|>orl develops ami discusses s t a t i 1 ic al methods 
lor processing nuclear attack information that have been mentioned in the 
preceding chapter. These methods are developed primarily in terms of a 
standardized report model that utilized time o! burst, reporting station 
location, azimuth, and range between reporting station and detonation. 

The model develops estimates lor the number of weapuns used, 1 he locations, 
times oi burst, and heights ol burst o1 the weapons, the weapon yields, and 
tiie targets against which the weapons were aimed. 

Knowledge of ground zeros and times ol burst is a prerequisite for 
ei i latent allot ut ion ol recovery i (.‘sources and is also important in pre¬ 
dicting fallout paths. Out a on ground zeros (together with data on yields 

and fission ratios) can also lie used with damage assessment techniques to 
provide civilian and military planners with estimates of surviving national 
resources. The estimates ol intended targets can be used to predict the 
attacker's choice of targets. Such prod let ions might enable evacuation or 
movement to shelter to be earned out before a location is actually attacked. 

Although manual jYUDF.T and sensor h'UDET reports may be the principal 
information sources lor this statistical process, the method is presented in 
a general way that is not dependent on these systems, Thus, major changes 
in reporting systems will not require corresponding changes in the statis¬ 
tical processing. This independence is possible because the different 
reporting systems provide similar info motion even though the formats and 
methods are very different. 



* This chapter is based on the working paper Statistical Methodology for 
Nuclear Attack Iniormnt ion Processing," TN-OAI’-h’E, by Charles R. Thompson, 
Stanlord Research Institute, Menlo Park, Cuiii . , August 19(37. 
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The Report Model 


The standard report model consists of data that can be obtained or 
Inferred from a report about iui attack. In some cases, the data may be 
available directly from a report by a station, but this may not necessarily 
be the case. 

The items of data in the standard model are: 

t = azimuth of burst from reporting station 
R -. range of burst from reporting station 
t - time of burst 

k - identifier of the reporting system used 
g = geographical coordinates of the repotting stution 
y = weapon yield 
h =. height of burst 

Azimuth is given in degrees clockwise from true north. Range is in 
nautical miles. Time of burst is given in hours after H-hour (start of 
attack). Item k identifies the reporting system by a code such as 1 for 
t._. .-1 rCJQuT auu 2 foi sensor NUUfcT, 1 he coordinates, g, of the reporting 
stations are given in terms of the Universal Transverse Mercator system 
in hundreds of meters. 

The model describes the statistical properties of the report items by 
means of probability density functions and confidence regions. The range, 
azimuth, and burst time estimates will he assumed to be independent, nor¬ 
mally distributed random variables with means at the true values and known 
standard deviations dependent on the reporting system. 

The standard errors arc in hours, degrees, or miles; they may be the 
same for all stations in a reporting system, or may be derived from a var¬ 
iance given in the original report from a station. 

Confidence Intervals 


The confidence interval lor the burst time of a weapon reported by a 
station, is obtained from tables of the normal distribution. A number, b, 
is obtained from the tables so that 
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1 - Q 


w he re 

Z is t he standardized nonn.d variable 

Q is the probability ol a variable exceeding the confidence 

Tile confidence interval lor the true burst line, t ol the 
burst is 

t - be s t ' "■ i ! | >T 


where 


o is the standard deviation ol the station report < f burst 
t is the reported time oi burst. 

11 the colli idem c level were 'Jit parent, then 

t - 1.0 1~ o t ' a t -> 1 (i1;» 

would be the confidence interval. I < > r a large group ol attack l 
inputs, y(j percent of the intervals const rue ted in this manner u 
the time of burst. 

Confidence Regions 


The confidence region for 1 he station-reported location of 
zero (fi, r) can alse be developed. One region i an lx; del mod by 
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0, - standard deviation ol azimuth 
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This region Is a sector of a circle about the 1ine-of-sight determined by 
the reported azimuth. The constants, a and b, tire obtained from tables t 1 
the normal d 1 s t r 1 bu t loti, so that 

r I r 1 

P r -j - a Z s. a 1 :< P f . | - b s. 7. s b| = 1 -cr 

If the confidence level Is TO percent, then a - b = 1.960 will proit.-.e i 
suitable confidence regio.. for (0, r), the true position of the ground 

zero relative to the reporting station. For a large group of reports, 

90 percent of the Intervals constructed In this way will cover the true 
ground zero. 

Figure 14 lllustrutes 90 percent confidence regions for the location 
of an actual ground zero by two different reporting stations. The con¬ 
fidence regions for the two station reports overlap, indicating that a 
single burst may have elicited both reports. Note that a confidence re¬ 
gion is associated with each repo-t. 

i ne commence regions may be approximaieu by rectangles because the 
estimated ranges, R, of the bursts will be large with respect to the range 
standard deviation, and small variances in angles are to be expected. The 
rectangular representation will allow manual processing of station reports 
by graphical techniques, as a back-up for a computer processing system. 


Associating Standardized Station Reports with Burst: 


The standardized station reports are the basic data for determining 
weapon characteristics. Each report gives data for a particular burst. 
Moreover, the same burst may have been reported by many stations. The 
more stations reporting a given burst, the more data available to develop 
information on weapon characteristics. In order to develop this informa¬ 
tion, the reports must be separated into groups corresponding to separate 
bursts. Because of the statistical nature ol the separation, it will be 
impossible to tell exactly which burst a station reported, thus causing 
some reports to be associated with more thun one burst. 
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The ron f idence region;- nmi time- o( burst associated with the station 
reports can serve as elass 1 fieations of reports into groups pertaining to 
u single burst or flusters oi bursts. The bus 1 1 ideas used m this classi¬ 
fication are: (i) the same t'lml Idence region tor several reports indicate 

a single burst for ail; (2) diilerent confidence regions lor two or more 
reports indicates diiterent hursts; and ( 3 ) overlapping oi confidence re¬ 
gions indicates that the reports may be for the same burst. The time of 
burst intervals may then be u-rd as a choc's oil the categories set up by 
confidence regions. For example, ll confidence regions overlap for two 
reports, then overlapping of the confidence intervals for their 'must time 
would indicate that 1 ho two reports pertain to a given single burnt, and 
they would be placed in that group. Figure 11 shows the overlapping ease. 

Detonation Hecurd Grouping 

The simple rule stated above is complicated by the "chaining' of con¬ 
fidence regions. figure lf> gives an illustration oi thi- dui.i. to el i« , 
reports 1 and 3 are .separate, but they are chained together lay report 2 . 
There may be a single burst in region if or bursts it, region; 1 and 3 . 

Since it is undesirable to ivcre.s l l myi t e tin number of bursts, a single 
group- containing records 1, if, nr.d it is made, subject to separation by 
times of burst . Similarly, records >1 to 8 also form a single group. 


Burst firm. 


The records grouped In posit lull may lx separated further on the basis 
of times of burst. Ill some- eases, the 1 1 mo.- may differ widely and hence 
allow a simple discrimination. Fur example, ll report b of Figure 15 has 
a burst time of b hour; later than the other records m t tic group, and the 
standard deviation is 1 hour, then record b must represent a distinct 
burst. After the removal oi record b from the group, n cheek indicates 
that two position groups are now present. In the absence of further time 
discrimination, three burst groups are formed from the original group of 
five records. They are ((>) , (1, 5), and ( 7 , 8); e.nh represents a dis¬ 

tinct burst. 
























Ill many cases, the times of burst intervals will overlap, but dis¬ 
crimination can still be performed by eonstruiting frequency plots of the 
time data. The modes of these plots can be used ns estimated burst times, 
and the records that have confidence intervals containing these estimates 
are grouped together. 

The plots are obtained by sorting the upper and lower bounds of the 
confidence intervals by time and assigning frequency values from the lowest 
time to the highest. If the sorted values are numbered by i - 1, 2, ..., n, 
and if f^ is me frequency of the ith value, then 

1 \ = 1 

fj=fj+lifiisa lower bound 
= f x - 1 if i is an upper bound 

It s 'ouirt be clear that f n - 0 because there are equal numbers of 
upper and lower bounds. 

An example will clarify this process and also illustrate the reduction 
of position groups by time of burst. Table 3 gives the estimated burst 
times from station reports 1 to 8 (Figure 15). If a 90 percent confidence 


Table 3 


REPORTED TIME ft OF BURST FOR RECORDS OF FIGUPE 15 


Record 

Time of Burst 
(hours after H hour) 

Standard Deviation 
(hours) 

1 

3 . 0 

1 . 2 

2 

3.5 

1 2 

3 

10. 1 

1 . 2 

4 

8.2 

1.2 

5 

1 (). 1 

1.2 

6 

13.1 

1 . 2 

7 

0.3 

1.2 

8 

7. 1 

1 . 2 


00 






level is required, then the confidence interval about each reported time i s 
t. 2.0 hours. The sorting ol record group (1, 2, 2,) will produce (l, 1.5, 

5, 5.5, 8.1, 12.1); the frc<~ucncy plot iron; this sorting is given in 
figure 16. Note that the frequency goes lu zero before the records are 
completely plotted. This indicates that the records should be separated 
into two groups: (1, 2) and (2). 

Figure 16 also shows a plot of records 4 to 8. This plot does not 
provide the high confidence separation that the plot of records 1 to 3 
provides, but it does have two modes. The modus at 12 and H hours indi¬ 
cate that two groupings should tie formed. They are (4, 5, 7, 8) and (5, 6). 
The first group can be further subdivided on n location basis into (-1, 5) 
anil (7, 8). This illustration also shows the use* of a single record, 5, in 
two different burst groups. 

Hurst Time Pi s<. r l rn a.; t : on wi th More Than Ore Reporting System 

The records represented in the above examples (Table 3 and Figures 15 
and 16) are a Mmplil ied case since all records ; me from the same* reporting 
system (standard deviation 1.2 hours). In practice, at least, tvo reporting 
systems (two cli1ferent standard deviations) will be in use; any set of 
records will be a mixture of records 1rom various systems. Figure 17 illus¬ 
trates this general situation tor two reporting systems. 

The rules for separating those i‘-cords into groups are similar to those 
for the simple case. The smaller cotif idence areas are grouped w-ithout ref¬ 
erence to the larger areas, Any larger areas overlapping the small area 
groups are added to these groups. Alter this operation, any larger areas 
not yet processed are grouped, and these larger ureas that overlap these 
groups are placed therein. For Figure 17, the following groups are derived: 

Croup Records 

1 9, 1, 2 

2 10, •}, 5 

3 11, 12, Id, 8 

4 d, 2 , 4 

5 6, 7, 8 
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Some of the records appear in more than one group because of the uncertainty 
of an exact placement or grouping, 

The position groups are tested for time of burst discrimination by the 
method of frequency distributions. The confidence intervals are of differ- 
ent sizes when different standard deviations are involved. The frequency 
function is modified to account for this by weighting the f^ inversely to 
the standard deviation. The modified rules are: 

fl = lo l 

*i - + 1/tfi if i is a lower bound 

- f x - 1/Oj if i is an upper bound 

The bounds have been arranged in order of reported times and i-numbered 
by i = 1, 2.n . 

Estimating Actual Ground 2ero 


Each group is considered to represent a single ground zero. The co¬ 
ordinates of the ground zero are unknown, but may be estimated by statisti¬ 
cal methods using the azimuth data from the standard report. Range data are 
not used because they are generally less accurate than azimuth data. If a 
burst group contains a single report, then the range and azimuth given in 
the report determine the estimated ground zero. 

A method for determining the position of an object from several 
azimuth measurements has been developed and presented.* This method can be 
used to determine an estimated ground zero which minimizes the sum 


• -1 ■ 


1=1 


* John, Floyd I., "Statistical Problems in Position-Fixing," WP-63-5, 
Office of Research Analysis, Office of Aerospace Research, Holloman 
Air Force Base, New Mexico, May 1963. 
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where 


= reported azimuth ol the l' ^ report 

= azimuth ol the estimated ground zero Iron the 
reporting station 

Oj = standard deviation of reported azimuth for the i 1 * 1 
reporting station 

n =. number of reports In the report group 

Figure 18 illustrates the general idea behind the method. 

In order to apply the method, a rectangular coordinate system must be 
established. Let (X^, Yj) be the coordinate of the reporting station 

in this coordinate system. The lines oi sight from the i*^ reporting 
station to the burst are written in terms oi this system as 

Y cos - X sin 9j = P 1 


where 


5 


i 



(X, V) = coordinates of an arbitrary point on the 
line of sight 


P = Y, cos Gi - X i sir. 9. 

The appropriate transformation for establishing this system is a 
gnomonic centered near the reported ground zeros. The earth coordinates 
of the station and reported ground zeros are moved by this transformation 
to a plane tangent to the earth at the center of the transformation. The 
projection is along a radius of the earth. The important feature of the 
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E s tima11ng the Target 


At 1 lie time of launch or release, every weapon has an intended target. 
The standard reports van be used t> provide information regarding the 
target to be associated with eaeh burst group. 

I.et the important resources that would be considered as targets be 
ldcitij- led by T, , T... , . . ■ lj„. i'ui e.u.u oi Luese target v uiiti lua ie>, i lie 

rdi;iz.• e must be hi.own in order to calculate azimuths and ranges with 
respect to the reporting stations. Targets with large areas for which no 
one set ol coordinates is satisfactory must be considered separate!!. 

The burst groups and target candidates are related by the statistic 


S 




u here 

n - number ol station reports in group 

;j, R. = station reported azimuth and range to estimated 
ground zero 

‘ , r - calculated azimuth and range to target j from 

i.l i.l 

reporting station 

r , q^ = azimutn and range standard deviations 

under the hypothesis that the intended target was T , the statistic S- 

J J 

is distributed as a chi-square variable with 2 degrees of freedom because 
the parameters are estimated from sample values. The hypothesis can be 
accepted or rejected by reference to the distribution tables of the chi- 
square variable. 

In some cut' , all target candidates will be rejected—considered as 
not being associated with a burst group. This leads to three possibilities: 
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1. Tho i''.' was tut burs I for the burst group 

2. The target candidate list is incomplete 

3. The weapon missed its target by a large distance 

Heuristics may lie developed to resolve t lu "e questions i r. operating 
sysferns . 

A further possibility is that 1 lie target was a resource area of 
large size, such as ail urban area or oil 1 leld. The reported ground 
zeros of the record group can be tested against t he area targets. An 
area may be accepted as a target lor a burst group i* a high fraction ol 
the reput i s in th>' £*’<»!•;• 1 -. 1 ' eatr a ground zero in the area. 

It may also occur that two target' are accepted lor a single burst 
group. This may mean that there is more than one burst represented by 
the group, b.it it may also mean that the data are not good enough to 
discriminate targets in all cases. A useful heuristic that can be used 
m this case is to weight the two targets by strategic importance and 
resolve ambiguity by selecting the more important target. Many variants 
ot this weighting system can be developed. One such development has 
been presented by Launno’ as a likelihood function, and provision for 
such functions is made in the subsys untt ol figure tv Initial Target 
Select ion . 

Until heuristics or more advanced statistical models are develoj>ed, 

however. a rule t ha t uses the stati tics j> ( j = 1,2, . . . , tr.) is used to 

j 

estimate the targets oi an attack. The rule is to select every resource, 

T^, lor which S^ is significant at the 10 i>ei'cent confidence level. The 

selected T are the estimated targets ol the attack, 
i 


* Launno, Richard, David Goodrich, and Donald Doane, "Statistical Methods 
of Estimating Damage from Nuclear Attacks on CONUS <U)," Stanford Re¬ 
search Institute, Menlo Park, California, December 1962 (SRI 2-1506, 
SECRET). 
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STAT1 STIt'AI. Sl'IiVIVAI. I- Si 1 MAT I NC CONI’KIT FI.O'.V CMAItT.S < OKTA 1 U,l>) 

Section VII ol t his r> port , the various .subsystems * ompr i s i ng the 
s' • .si leal survival esi 1 mat nig. c uni ept are described, and I low (.[tails 
fire iresented which outline Use steps m each subsystem activity. Tins 
appendix contains more detailed 1 low chart.' lor some of the subsystems. 
These more detailed 1 low charts were devised dus lug development ol the 
basic methodologies presented by this repori ; as such, they are cons id- 
eied to be over and above tin basic objective ol this study. Tl.t.y repre¬ 
sent it tirst step in i ntp Icn.eti t i ng the l>a s 11 me t bodo I og les tor computer ir-c 1 , 
and they are documented here so that the eI Io> t expended will not be lost, 
(inures A-l, A-2, 1-3, and A-1 present 1 low charts lor the following 

He t 1\3 tles ; 

1. Data Input and let'll n it ion 

2. Pel ini lion ol the Initial I at got Systen, 

3. Pelinilion ol the iinul Taigct System 

1. t’arainet in Survival K~ 1 irut t i ng 

1 liese 1 low charts do not iollow the same sequence as the subsystem 
a. l i \ i t i e s tie s e i 1 bed 1 ii Se ( t 1 1 u i ill. 1 oi (.- .cum p re , i lie <ict i l i t i uutl i neii 
m figure .1-1, Data Input and Vi'i'i I H at lull, includes most ol tile activi¬ 
ties described in Section VII under the headings: data input and verili- 
eatioii, direct damage report processing, requested status report proc¬ 
essing, and detonation icimd maintenance. Also, tItc level of detail is 
not ci n i f o l in throughout the 1 low charts. This rellects the fact that at 
present, certain activities within the proposed statistical survival es¬ 
timating concept have received more study, and are more sjieel t ieal 1\ de¬ 
fined, than other activities. 

The concept presented in Figure A--1, Parametric Survival Ks t ima t i ng, 
represents a somewhat dil'leri-nt approach to parametric survival estimating 
Iroiii ilia I given in Section I'll. In general, it is considered by the 
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authors to be aoru cumbersome anil probably no more accurate than the 
paramo trl r survival estimating procedure of Sort Ion VII. It is included 
In' re only for possible further evaluation when the busic methodologies ol 
this study are Implemented. 

In spite of tlulr 1ncompleteness, the following llow rhnrts should be 
uselul uecnuse they represent a lirst attempt at defining the logic ami 
interrelationships of certain essential activities of the statistical sur¬ 
vival estimating concept In a more detailed manner. With additional work, 
they could be extended to provide the basis for a computerized system imple¬ 
menting the proposed s' 1 '-' t ■"-* 1 estimating concept, 
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Arr-eudts i: 

UNO. hi hi idNA' ion i'.toy. 

inUKO DAMAliK lU.PUKTS 


ha r I ill se« t ions "■ 11;».v i< jii'it have stated liul dcton.it ion i icuril.i 

may bn ini nrred from direct damage reports and I i-uiu sUlreports. The 
prim-i|>;!l (tala that must hi interred tin hide weapon yield, him:lit tj I hurst , 
ami ground «-ni, It is hoped ih.it mri-. t o.image iup'i't- and st.it u- reports 
will contain time ot Inn st ini ornttii ion. Jf not, however. time ol burst 
must also be interred. and in this ease t be time that the it poll was made 
would probably become t lie principal basis lor itileii inn time ol burst. 

Had in)onicai fallout inlormattun nnpiu also be uselul for Ibis purpose. 

Various techniipies and me thuds aii a\.illable tor inferring t lie weapon 
yield, height ol burst, ami ground /„ci i data. A method lor estimating 
these essential data from radiological iallmit reports has been devised 
and piesented b\ Uuthun.* Other : ...llile method' nu lude pattern rccog- 
n 11 ion. hypothesis testing, statistic a 1 applications, and heuristic deei- 
s t on t ii 1' -, 

A s\stem based mainly on pattern rccognition has been suggested by 
11 ie 01 1 i ce of Cl vi 1 liei cjist, V. i tb. 1 1,iv aystent. cmergcaey ojx rat i n g center:, 
will tins; receive estimates ol ope rat itig situations from subordinate op¬ 
erating areas and plot this information on a map to establish the general 
situation within that region. The basic operating situations corresponding 
to the suggested map codes (1 through 9) are shown in Figure ii — 1. 


* Hot him, li. C. r Uses oi Radiological Fallout Information in Oi* rational 
Post a t t 8ck I‘a mage Assessment Systems (li), st aidm-il Research Institute, 
Menlo Puk, Calilornia, January 19t>f> (CONFIUKNT1A1.) . 
t Concept ol Operations under Nuclear attack," Ulllce oi O u I lieleiise, 
Washington, U.C., June 2b, 1907, (Working l.'ralt) 




















Figure 8-i 

NINE BASIC OPERATING SITUATIONS 
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Source: "Concept of Operations under Nuclear Attack," Office of Civil Defense, 
Washington, D.C., June 26, 1967 (Working Draft) 
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Figu '2 A-l 

DATA INPUT AND VERIFICATION 
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Figure A-2 

DEFINITION OF THE INITIAL TARGET SYSTEM 
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Figure A-3 

DEFINITION OF THE FINAL TARGET SYSTEM 
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Figure A-4 

PARAMETRIC SURVIVAL ESTIMATING 
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It is suggested that early reports from IAIFIKK amt adjacent NEGFIKE 
/.oner^ in the vicinity ot a nuclear explosion could then be used to estab¬ 
lish the location and approximate yiul1 ol the detonation by measurement 
of 1 he diameter of the LOFIHK area and comparison with a curve or table 
showing the anticipated diameter as a inaction ol yield. An alternative 
procedure would be to apply templates ol predicted 1.0F1KE areas for weapon 
yields to a map showing operating area situation reports to estimate loca¬ 
tion and yield ol detonations. The later procedure is illustrate., in Fig¬ 
ure H-H. In this illustration, the I.OFIUK region extends to approximately 
the 1 psi range, which is in agreement with generally accepted estimates 
I or si at t e red 1 lres. 

A statistical approach to the estimation ot weapon yield is also pos¬ 
sible. One possibility is to make u>c ol small sample theory techniques 
such a-, use ol the noncentral "student's" t distribution, which tests the 
hypothesis that the mean perimeter ojierating area radius is equal to a 
speciiitd value, Jn tills distribution, t -tuti-ti- t 1 - du lined as: 


V N - 1 


N number ot sample point-, or peri meter o|>eruting areas. 

n = hypothesized mean radius be tween ground zero and centers of 
Ijci'iraetcr operating area-, i .«•, . foi 1 MT, a MT, 10 MT, or 

20 Nrr. 

X mean ol sample radii between ground zero and centers ol perimeter 
operating areas. 

s - standard deviation ol sample radii. 


This system might be used in such a way that u lor the weapon yield is 
tested 1 or various hypothesized values ol |ierhaps 1, a, 10, and 20 mega¬ 
tons. The use of a statistical approach has the important advantage that 
coniidenee intervals may be set which supply the system operators with in¬ 
formation as to how much coniidenee might be placed upon a given weapon 
yield seleclion. 













Figure B - 2 
NUDET ESTIMATE 



Source: "Concept of Operations under Nuclear Attack" 
Office ot Civil Defense, Washington, D.C. 
June 26, 1967 (Working Draft) 
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The statistical approach might also be used in estimating ground zero; 



f 

for example, a circular area of weapon effects may be derived through pat- 

H 



tern recognition for two or more of the weapon effects such as observable 



i; 

blast damage and scattered fires. A comparative analysis of these two 



i 

circles for the purpose ol estimating ground zero is possible that would 

- 



have a more quantitative basis than simple pattern recognition. In actual 
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i 

practice, however, It mav well be that simple pattern recognition Is ude- 

_£ 
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quale for estimating ground zero. 
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The above discussion has illustrated that in many oases, a combina- 
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tion or pattern recognition, heuristic decision rules, and statistical 

| - -- 



analytical techniques can be used t u estimate essential NUIJKT data 1 ront 



I 

extremely elementary direct damage reports. Yet because ol the obvious 

i - - -^z 
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uncertainties involved, sensor and manual NUDKT reports should l>e used 



i 

when available. An obvious difficulty with the inference of detonation 


t 

- 

reports 1 rout direct damage reports is that direct damage reports probably 
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jr 

will r.ot have been received from all operating areas, particularly in the 


r 


ear 1 \ time periods. When insufficient direct damage reports are available 



i i 

for estimation purposes or when additional direct damage reports are de¬ 
sired. status reports should be requested for the appropriate o;>erating 
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A STATISTICAL INFORMATION ST5TEM FOR ESTIMATING 
THE MAGNITUDE AND SCOPE OF NUCLEAR ATTACKS 
by Robert M. Rodden 
Stanford Research Institute 
February 1968 

OCD Work Unit 4615A 
DETACHABLE SUMMARY 


This research develops and presents a statistical concept for making 
national survival estimates in the transattack and early postattack 
periods, Statistical methodologies are used to process reports ot detona¬ 
tions, and the principle of statistical inference Is used to estimate the 
magnitude and scope of the attack. Existing systems for the preparation 
of national survival estimates in the trarsattack and early postattack 
periods are based largely on preattack planning methods or on the tabula¬ 
tion ci direct damage reports as they are received. Adequate and accurate 
dat; of the kind needed to produce survival estimates with these techniques 
will not be available in the transattack and early postattack periods. 

The lour principal activities of the concept developed by this re¬ 
search are (P real time input processing and data preparation, (2) sta¬ 
tistical conversion of data to information, (3) inference of attack size 
and target system, (4) estimation of effects on resources and population. 
Activitie. 1 and 4 above may have much in common with certain existing 
survival estimating systems. The main distinctions between this and other 
survival estimating systems are in activities 2 and 3. The concept devel¬ 
oped by this research uses statistical inference and heuristic decision 
rules to estimate the magnitude of the total attack. The methodology thus 
devised can operate with incomplete and inaccurate information. Targets 
are interred by statistical inference alter appropriate grouping of tar¬ 
get candidates into target categories and vulnerability groups. 

The methodology Includes: the receiving and processing of burst in¬ 
formation and direct damage reports for use in making survival estimates; 
basic concepts for organizing inputs from various sources into appropriate 
data files; and provision for supplying measures oi reporting performance 












to ihe system operator. Detonation reports and other data are converted 
statistically to information that is more directly applicable to survival 
estimating. Estimated actual ground zeros and targets for the initial 
target list are determined statistically. Areas whore data are incomplete 
or possibly inaccurate are identified, and status reports are requested to 
correct these deficiencies. An initial target system is developed that 
yields information on the target categories hit and on the severity of the 
a t tac k. 

The next main step is inference of the full attack size and target 
-system. This is accomplished with the aid of the initial target list, a 
list of target candidates, and a knowledge of target categories and vul¬ 
nerability groupings. At the completion of this step, an augmented target 
list that in general will be expanded considerably from the initial target 
list will have been developed through statistical inference, Again, 
status reports will be requested to cover areas of doubtful or inadequate 
i n formation. 

The final step of the concept is to estimate ellects on resource^ and 
population by means of the synthesized target system. It is shown how a 
target system may be coupled with several survival estimating techniques 
to produce survival estimates. The survival estimating technique selected 
may well depend upon the requirements ul the system operator at a given 
point in time. A summary ol the system and its principal activities is 
given in Table 1. 

Section VIII presents details of statistical methodology for esti¬ 
mating actual ground zeros and targets relevant to a group of detonation 
reports. A basic tool is the use of coni idence regions and coniide.ice 
intervals to determine the actual number of bursts and estimate actual 
ground zeros associated with a particular group of detonation reports. 

The selection ol target candidates for the initial target list is also 
accomplished statistically, using a least-squares methodology. Weapon 
yields and heights of burst lor a given group of defoliation reports are 
derived mathematically. 

Appendix A provides mere detailed flow charts for some parts of the 
conceptual system. The basil purpose ol these flow charts is to document 
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ther ideas that have been developed in the course of this research but 
that are not appropriate to the generalized flow charts of the main re¬ 
port. The ilow charts of Appendix A will help provide the basis for de¬ 
veloping the system to a point where computer programs for the system can 
be written. However, considerable additional work will be needed before 
these flow charts are adequate for that purpose. 

The feasibility of inferring detonation records from direct damage 
reports is investigated briefly in Appendi.\ B, which also gives a frame¬ 
work for such inference, using pattern recognition and hypothesis testing. 

This research effort has develojied methodology for estimating attack 
size in the transattack and early postattack periods, and has shown how 
this methodology can be combined with existing damage assessment methods 
to produce national survival estimates. More advanced survival estimating 
techniques, that would integrate attack size estimating and damage assess- 
nen' more directly, are briellv explored. Survival estimating systems 
recommended for lull development and implementation by the National Civil 
Defense Computer Facility are identified. 
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